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общей для всех методов тестовой задачи исследования много- 
экстремальной функции Растригина. Показано, что все ука- 
занные методы вполне пригодны для решения многоэкстре- 
мальных задач. Хотя в каждом из эвристических алгоритмов 
приходится использовать собственные специфические подхо- 
ды к решению задачи обнаружения и идентификации локаль- 
ных экстремумов, их объединяет необходимость осуществле- 
ния кластеризации данных. Каждый метод может обеспечить 
любую заданную точность решения экстремальной задачи и 
использует приемлемый ресурс времени. 
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роящихся частиц, муравьиный алгоритм. 


Введение. Многие современные проблемы науки, техники, экономики, военного дела и пр. связаны с решением задач 
поиска оптимальных характеристик объектов проектирования: конструкций, технологий, режимов и условий работы, 
динамических и статических состояний и т. д. Иными словами, разработчикам приходится решать задачи поисковой 
оптимизации (ПО) [1-3]. Характерно, что большинство известных на сегодня методов ПО разработано и эффективно 
используется для нахождения одного оптимума — чаще всего, глобального [3, 4]. Однако многие задачи планирова- 
ния, сложные технологические комплексы, транспортные задачи и другие объекты оптимизации (особенно дискрет- 
ной природы) характеризуются многоэкстремальностью (МЭ) [4-11]. Столь существенное отличительное свойство 
требует специфических методов решения таких задач. Вряд ли эти методы целесообразно искать в классе детермини- 
рованных методов ПО. Они слишком чувствительны к знакопеременности и разрывности функций отклика в конти- 
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нуальных факторных пространствах, а также описываются МР-полными алгоритмами в дискретных факторных про- 
странствах. Для решения большинства реальных оптимизационных задач все чаще стремятся применять методы, по- 
лучившие название «эвристические». Эти методы наиболее перспективны и для решения МЭ задач [5—11]. 
Постановка задачи. В связи со сказанным во введении представляется актуальным и перспективным исследование 
наиболее распространенных эвристических методов ПО в среде какой-либо характерной МЭ задачи. Цель работы: 
установить возможность нахождения нескольких или всех экстремумов рассматриваемой области каждым из выбран- 
ных методов. Кроме того, необходимо оценить точность определения значений экстремумов и их координат. Поэтому 
первым этапом исследований является выбор тестовой МЭ функции, которая может послужить общей для всех мето- 
дов средой решения МЭ задач. 
Выбор многоэкстремальной тестовой функции и предварительный анализ ее свойств. Самыми распространен- 
ными и эффективными тестовыми функциями, используемыми для отладки и исследования методов ПО, являются 
функции Розенброка, Химмельблау и Растригина. Из них наиболее многоэкстремальная — функция Растригина (ФР). 
Эта невыпуклая универсальная функция предложена в 1974 году Леонардом Растригиным [12]. Уравнение функции М 
аргументов имеет вид: 

Ш 

1%) = Ахп+ У[2 — Ахсоз(2хлхх)] (0 

1=1 

где = 10. 


Глобальный минимум этой функции — Г (0)=0. 


тт 
Нахождение глобального минимума данной функции представляет собой достаточно трудную задачу из-за 
большой области поиска и количества локальных минимумов, на которых алгоритм может зацикливаться. Выделение 
и оценка локальных экстремумов является еще более сложной и неординарной задачей. 
Исследование ФР алгоритмом роящихся частиц. Сущность и основания для использования в задачах ПО метода 
роящихся частиц (МРЧ) хорошо известны [13—17]. Классический алгоритм МРЧ имитирует реальную групповую по- 
веденческую модель насекомых, птиц, рыб, многих простейших и т. п. Однако МЭ объекта оптимизации требует не- 
которых специфических свойств алгоритма. Поэтому канонический вариант МРЧ был значительно переработан авто- 
рами [18—20] и другими учениками Р. А. Нейдорфа [7-11]. В частности, разработана его модификация для решения 
МЭ задач. 
Модель механического движения (ММД) частиц [20] в МЭ МРЧ была существенно доработана и уточнена: 


Хи= Ха-му+ Рам ХДЕ, (2) 

: = Рам ы: Ам ХДЕ, (3) 

А:= Ар:+ А (4) 

где Хм); — предыдущее положение 1-й частицы; Х’,,— ее текущее положение; 7, — скорость 1-й частицы в теку- 
щий момент времени; Ре: ^№): — предыдущая скорость частицы; Яр №): — предыдущее ускорение частицы в преды- 


дущий момент времени; ДЕ — интервал интегрирования по времени; Ан ускорение под действием биологически 


обусловленных сил притяжения частиц; А,.; — замедление под действием сил трения. 


т 
Для улучшения ПОИСКОВЫХ свойств было введено стохастическое размытие параметров: 
(в) =Ах(1+2х&(тта()-0,5)), (5) 
где ^5( 5 ) —Щ_ значение флуктуирующего параметра на такте вычисления, & — относительное отклонение искажаемо- 
го параметра от номинального значения; ги4(1) — случайное число в диапазоне [0, 1]. 


Предложены и проверены механизмы естественной кластеризации частиц, порождаемой и параметрически 
настраиваемой свойствами алгоритма: 
— градиентный, основанный на чувствительности частиц к смене знака скорости изменения оптимизируемого крите- 
рия исследуемого объекта [7-9]; 
— потенциальный, основанный на введении в ММД сил притяжения ко всем обнаруженным при роении и сканирова- 
нии пространства поиска локальным экстремумам: 


Ар = Ар. + Ар + р . (6) 
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-с . ЗЕ р 
Здесь А»; — притяжение частицы к глобальному экстремуму; А’; — притяжение частицы к локальному экстремуму; 


я — притяжение частицы к центру ближайшего кластера. 


Алгоритм МЭ МРЧ показал хорошую избирательность при локализации областей экстремумов, однако меха- 
низмы кластеризации и локализации кластеров потребовали существенной структурной и параметрической доработ- 
ки. В настоящем исследовании введены и проверены следующие модификации: 

— механизм отсева по определенным критериям «плохих» кластеров (являющихся худшими на заданном числе ите- 
раций); 
— механизм объединения схожих кластеров на каждом шагу; 
— настройка условного притяжения к центру ближайшего кластера; 
— механизм локализации областей кластеров и поиска в них параметров локальных экстремумов. 
Эффективность модификаций проверена ФР в диапазоне координат (х, у) в [-1.5,1.5]. В этой области ФР 


имеет 9 локальных минимумов, один из которых является глобальным. 
На рис. 1 отображен процесс локализации областей экстремумов и создания соответствующих кластеров. 





а) 6) 
Рис. 1. Локализация областей экстремумов: 1-я итерация (а); 50-я итерация (6) 

Как видно из рис. 1, частицы изначально притягивались к образовавшемуся кластеру, который располагался в 
области глобального экстремума. Затем периферийные частицы находили локальные экстремумы, притягивались к 
ним и собрались в кластеры. В суженных областях кластеров применение алгоритма МЭ МРЧ (на случай наличия ме- 
нее выделенных и значимых экстремумов) повторяется. Этот процесс итеративно повторяется до получения требуе- 
мой точности оценки параметров всех локальных и глобального экстремумов. 

На рис. 2 отображены фрагменты итеративной идентификации локального экстремума, который располагает- 


ся в точке [-1, 1]. 





Рис. 2. Локальная идентификация одного из локальных экстремумов: 15-я итерация (а); 33-я итерация (6); 50-я итерация (в) 
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В табл. | отображены результаты, полученные локализацией всех областей, приведены координаты х =х; и 
у =х», а также значения ФР, полученные по формуле (1). Как видно, получены удовлетворительные результаты с вре- 


менными затратами -—149 с. С увеличением числа итераций (и времени поиска) точность оценки повышается. 









































Таблица 1 
Результаты проведенного опыта 
Эталон Параметры оценки экстремума 
Координаты Значение 
х у Л(х, у) ы - У 
—1 1 2 —0,9957 0,9953 1,9901 
—1 0 1 —0,9949 0,0001 0,995 
—1 —1 2 —0,9951 —0,9947 1,9899 
0 1 1 —3,20х10° 0,9948 0,995 
0 0 0 9,8510 —6,49*10° 1,94*10°° 
0 —1 1 0,00017 —0,995 0,995 
1 1 2 0,9949 0,995 1,9899 
1 0 1 0,995 0,00011 0,995 
1 —1 2 0,9952 —0,9951 1,9899 























Для решения задачи МЭ МРЧ использован ПК с процессором АМО Р®|епот П Р960, оперативная память — 
6ОЪ. При этом для достижения точности 1Е-3 удалось ограничиться 50 итерациями и осуществить локализацию обла- 
стей экстремумов с затратами по времени — 149 с. При дополнительном исследовании каждой из областей потребо- 
валось в среднем — 2,5-3,5 мин на одну область и порядка 30 мин на все 9 областей. 

Таким образом, можно сделать вывод, что МЭ МРЧ является эффективным инструментом решения МЭ задач. 
Исследование ФР эволюционно-генетическим алгоритмом. Эволюционно-генетический алгоритм (ЭГА) является од- 
ним из наиболее востребованных инструментов решения оптимизационных задач [21—24]. Структура и основные опе- 
раторы ЭГА хорошо известны, а конкретные структурно-параметрические особенности зависят от предметного при- 
ложения. В частности, применение ЭГА для решения МЭ задач [25—28] требует дополнения классического ЭГА ин- 
струментами селектирования экстремумов — по типу (максимум, минимум), по величине, а также по координатам его 
положения в факторном пространстве. В данной работе развивается подход к селекции экстремумов, основанный на 
использовании одновыборочного критерия Стьюдента [27—29]. Сущность данного подхода может быть представлена 
следующим образом. При последовательном использовании ЭГА результаты, полученные в финишных поколениях, 
кластеризуются, т. е. разделяются по координатным группам. С этой целью для каждого результата проверяется 0- 
гипотеза. 

Математическая модель данного подхода кластеризации заключается в последовательном сравнении векторов 
со средним значением группы векторов 


У={Ди=у-т = Я, 
где и — количество векторов группы. 
Относительно заданной доверительной вероятности принимается решение о принадлежности вектора у суще- 
ствующему множеству. Для уточнения принадлежности необходимо рассчитать среднее значение длин векторов 
группы, с которой производится сравнение: 


Ду=У" Ду/п. (7) 


После необходимо вычислить стандарт отклонения длин векторов уже обозначившегося кластера: 


Зил = | (47 У2Ит-И, (8) 


Зин (9) 


По вычисленным величинам рассчитываются опытные значения одновыборочного критерия Стьюдента: 
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юНДУ-Ау| / буду}. (10) 

Если найденное опытное значение 1, не превосходит табличного [30] значения (. при п степеней свободы и 
выбранном в таблице уровне доверительной вероятности Р, то можно считать, что и принадлежит данной группе объ- 
ектов. 

Данный метод хорошо зарекомендовал себя при исследовании МЭ зависимостей [27-31]. 

С учетом этого подхода разработаны алгоритм и программное средство (ПС) на его основе. С помощью ПС 
ФР рассмотрена в том же диапазоне, что и в предыдущем разделе статьи при исследовании ФР алгоритмом роящихся 
частиц. 

Структура входных параметров ЭГА, которая была использована при исследовании ФР, включает в себя по- 
колений — 10, особей в каждом поколении — 1000, вероятность кроссинговера — 95 % и мутации — 30 %. Стоит 
отметить, что параметр области поиска используется тот же, что и в предыдущем разделе, а точность исследования в 
этой области равна 7 разрядам после запятой. В результате выделено 9 кластеров, минимумы которых можно соотне- 
сти с теми минимумами ФР, которые попали в исследуемую область. 

На рис. 3 отображены графики последовательного нахождения значений ФР и их различных координат (Х и 
У), а также соответствующие им значения целевой функции (Е (Х, У) = 2), которые, в свою очередь, упорядочены по 
убыванию (для кластеров, образованных вокруг минимумов со значениям (-1;1), (1—1), (1;1), (-1;—1)). 
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Рис. 3. Выделение кластеров в ходе эксперимента 


На рис. 4, а восемь периферийных кластеров характеризуют локальные минимумы, а центральный кластер со- 
держит результаты, приближающиеся к глобальному минимуму функции. Хорошо видно, что близкие (а в некоторых 
случаях равные) по значению функции цели имеют существенные отличия по координатным показателям (т. е. пара- 
метры функции цели, обеспечивающие близкие к минимуму значения, разнятся). Данный факт подтверждает много- 
экстремальность исследуемого объекта. 








а) 6) в) 


Рис. 4. Локализация областей экстремумов. Выделенные кластеры ФР (а) и формирование кластеров в локализованной области: 


100-е поколение (6); 110-е поколение (в) 
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Результаты исследования функции (как на глобальный, так и на локальные минимумы) отображены в табл. 2. 









































Таблица 2 
Найденные параметры функции цели первой итерации (по кластерам) 
Эталон Параметры оценки экстремума 
Координаты Значение 
№ х У Л(х, у) 
х ы. Л(х, у) 
1 0 0 0 0,00188 0,00015 0,00071 
2 —1 0 1 —0,98932 0,00073 1,00137 
3 0 —1 1 0,00824 —1,00043 1.01436 
4 1 0 1 0,99948 —0,01328 1,03398 
5 0 1 1 0,01403 1.00665 1,0611 
6 1 —1 2 0,99146 —0,99314 1,993 
7 1 1 2 0,99702 1.00467 2,00947 
8 —1 —1 2 —0,9897 —1,01397 2,06707 
9 —1 1 2 —1,00528 1.00077 2,01775 


























Как видно из табл. 2, оценки значений экстремумов и их координат не отличаются высокой точностью. В слу- 
чае, если найденные значения не удовлетворяют требуемой точности, реализуется повторное исследование миниму- 
мов функции в области кластеров. 

В качестве примера приводятся исследования областей кластеров 2, 4 и 8 методом локализации поисковой 
области [29] (табл. 3). Дальнейший поиск значений осуществляется использующим ЭГА методом, разработанным ав- 
торами [27, 28]. 

Наглядно кластер найденных ЭГА статистически близких значений в одной из областей можно наблюдать на 
рис. 4, 6. Наилучшая оценка экстремума выделена жирной точкой. Для повышения точности оценки вокруг этого суб- 
экстремума программой выделяется область дополнительного исследования, результаты анализа которой визуализи- 
рованы на рис. 4, в. 

Уточненные на второй итерации значения функций цели, а также соответствующих им координат приведены 
в табл. 3. 

















Таблица 3 
Найденные параметры значения функции цели второй итерации 
2-й кластер 4-й кластер 8-й кластер 
Значение Х —0,999996 0,999992 —0,999996 
Значение У 0,000005 —0,000008 —1,000005 
Значение ФР 0,999997 0,999989 1,999991 

















Исследование ФР муравьиным алгоритмом. Муравьиные алгоритмы (МА) — еще одна группа методов, которая ис- 
пользуется в различных оптимизационных задачах. Их отличительная особенность заключается в имитации ключевых 
черт поведения настоящих муравьев [32]. Чаще всего МА применяются к задачам минимизации пути на графах [33], 
однако данные алгоритмы показывают хорошие результаты и в других областях [34, 35]. В представленной работе 
рассмотрено применение классического МА для МЭ оптимизации эталонной двухпараметрической ФР [12]. 

В основе описываемого метода лежит классическая реализация МА, применяемая к задачам на графе [36, 37, 
38], но с некоторым дополнением. Рассматриваемое пространство делится на равные участки [39]. 

Как и в классическом МА, в данной модификации выделяются такие этапы, как «расстановка и инициализа- 
ция», «осуществление перехода» и «проверка условия останова». 

В качестве примера был рассмотрен фрагмент ФР для хе[-1.5,1.5], у[-1.5,1.5]. Он разделен на пхп 


участков, каждый из которых ассоциируется со значением функции в центре и наделяется некоторым количеством 
феромона. На участок помещается А муравьев. Так как все участки одинаковые, то размеры фрагмента можно вычис- 
лить по формуле: 


т = (Хтах- Хти//П. (11) 
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Таким образом, множество участков представлено матрицей А=(1,;; И ‚1: Муравей в участке Г,;,; оцени- 


вает все смежные с ним участки и рассчитывает вероятности перехода по формуле: 














р а Ур рьу, р > Лань ньд > 0%}; (12) 

ы У ьу дЛя Ь] БУнь > 0, 
где О*( ННьль» Ин р, ‚„0.Ви ) — функция количества феромонов в ячейках — “+1, /+1> 77; ПаулньЙ, = 
Уж, р )-У(жьнь а ) № — вес (виртуальная дистанция) между двумя участками; 0 — изменяемый коэф- 


ициент воздействия еромона; — изменяемый коэ ициент весового воздействия; { — номер итерации. 
) 


Структура функции О(тн нь Пн, нь И, „В, ) может быть различной [39]. 


Модель (12) дополняется моделью обновления феромона (13). На каждой итерации происходит прирост и ис- 
парение феромона. Таким образом, приращение количества феромона за шаг моделирования в участке 7, , рассчиты- 


вается по следующей формуле: 
х,(1+0)=(1-р)хт,(Е)+Алту (1) .(13) 
Здесь ре(ОЛ) — коэффициент испарения; 17 (1) — содержание феромона на участке / 17} Алу (1) — прираще- 


ние на каждой итерации, рассчитываемое по формуле: 
Ах; (Е о Кх Не (р) > а (кр) 14 
лу(. ) с ((хрУьР Лиль Уналны }) , ( ) 
рВ= 


где К — коэффициент прироста феромона. 
Испарение происходит в связи с тем, чтобы было рассмотрено все пространство возможных решений. 


При Л(х,у)и > Л(х,У),, переход на участок [,,,,;-1с участка /,, запрещен. Таким образом, условие 


останова считается выполненным, если все муравьи не имеют возможности двигаться. В результате через М итераций 
муравьи займут ячейки с наименьшим значением функции, локализовав минимумы. 

На основе описанного алгоритма и моделей (12, 13) разработано ПС, реализующее поиск локальных и гло- 
бальных экстремумов. В качестве примера на рис. 5 приводятся результаты поиска глобального и локальных миниму- 
мов ФР. Для решения выбраны границы области определения ФР, аналогичные принятым в предыдущих разделах. 
Выделенная область изначально делилась с шагом 0,25, и в каждом фрагменте помещалось по 2 муравья. Коэффици- 
енты а = 1, В= 0,5, р = 0,5, К =1ит= 1. На рис. 5 изображены отдельные этапы работы ПС. 











Рис. 5. Визуализация этапов работы ПС: инициализация (а); 3-я итерация (6); финальный результат (в) 


Уже после трех итераций наблюдается повышение концентрации агентов во всех областях минимумов функ- 
ции. Финальный результат получен за 6 итераций. На рис. 5, в видно, что агенты не только обнаруживают все мини- 
мумы, но и выделяют их значимость своим количеством. 

Рассмотрим пример применения алгоритма с большим разбиением исследуемой области. Данный алгоритм 
применяется повторно к локализованным участкам до тех пор, пока не будет получен результат с заданной точностью. 
На рис. 6 отображены полученные за 220 итераций результаты локализации и двух уточнений глобального экстрему- 
ма с разбиением 100х100 участков, который располагается в точке (0, 0). Такой результат получен в связи с исходны- 
ми данными и особенностями строения ФР. 
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а) 6) в) 
Рис. 6. Результаты локализации и уточнения при разбиении области 100х100: выделение всех локальных экстремумов (а); 
выделение наилучшего локального экстремума (6); уточнение глобального экстремума (в) 


В заключение стоит отметить, что при относительно малом количестве итераций МА дает прекрасные резуль- 
таты, что свидетельствует о его эффективности при МЭ оптимизации функций. 
Заключение. Для исследования эффективности решения задачи оптимизации при многоэкстремальности объекта бы- 
ли выбраны три эвристических алгоритма. Анализ результатов их использования показал, что данные методы дей- 
ственны и эффективны, хотя есть существенные различия между структурами процессов поиска и между формами 
оценки получаемых результатов. Эти свойства алгоритмов определяются в каждом конкретном случае специфически- 
ми подходами к решению задачи: обнаружение и идентификация локальных экстремумов, способы кластеризации и 
последующие операции уточнения результатов. Тем не менее, во всех случаях оказалось необходимым подвергнуть 
каждый из этих методов существенной модификации, связанной с необходимостью кластеризации данных. Исследо- 
вания показали, что каждый метод может обеспечить практически любую требуемую точность решения экстремаль- 
ной задачи (наилучший результат: МРЧ = 1,94Е-06; ЭГА = 5Е-06; МА = 3,75Е-07). Это, конечно, обусловливает значи- 
тельное увеличение вычислительных и временных ресурсов. Однако применительно к исследуемой тестовой функции 
все методы показали вполне приемлемую производительность (средние затраты времени на эксперимент: 
МРЧ & 300 с; ЭГА & 86 с; МА ®& 128 с). Таким образом, все три исследованных в данной работе метода актуальны и 
перспективны. Конкретный выбор алгоритмического инструмента для решения той или иной МЭ задачи зависит от 
опыта и от личных предпочтений исследователя, а также от специфических характеристик предметной области иссле- 
дования. 
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